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Abstract 

We investigate the combined effect of neutron and proton superfluidities on the 
cooling of neutron stars whose cores consist of nucleons and electrons. We consider 
singlet-state pairing of protons and triplet-state pairing of neutrons in the cores of 
neutron stars. The critical superfluid temperatures T c are assumed to depend on the 
density of matter. We study two types of neutron pairing with different components 
of the total angular momentum of Cooper pairs along the quantization axis (|mj| 
=0 or 2). Our calculations are compared with observations of thermal emission from 
isolated neutron stars. We show that the observations can be interpreted by using two 
classes of superfluidity models: (i) strong proton superfluidity with a maximum critical 
temperature in the stellar core T™ ax > 4 x 10 9 K and weak neutron superfluidity of any 
type (T c max < 2 x 10 8 K); (2) strong neutron superfluidity (pairing with |?nj|=0) and 
weak proton superfluidity. The two types of models reflect an approximate symmetry 
with respect to an interchange of the critical temperatures of neutron and proton 
pairing. 



* e-mail: gusakov@astro.ioffe.rssi.ru 
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1 Introduction 

At present, the properties of superdense matter in the cores of neutron stars are known 
poorly. For example, the fundamental problem of the equation of state for matter of supranu- 
clear density has not yet been solved. The existing calculations are model dependent and 
yield a variety of equations of state in the cores of neutron stars (Lattimer and Prakash 2001; 
Haensel 2003) with different compositions of this matter (nucleons, hyperons, pion or kaon 
condensates, quarks). The properties of nucleon superfluidity in the inner layers of neutron 
stars are also unclear. Calculated critical temperatures of nucleon superfluidity strongly 
depend on the nucleon-nucleon interaction model used and on the method of allowance for 
many-body effects (see, e.g., Lombardo and Schulze 2001). However, they can be studied by 
comparing the cooling theory of neutron stars with observations of thermal emission from 
isolated neutron stars. Here, we continue to simulate the cooling of superfluid neutron stars 
whose cores contain neutrons, protons, and electrons, with critical temperatures of nucleon 
superfluidity depending on the density of matter. 

We extend the class of cooling models that were proposed by Kaminker et al. (2001, 2002) 
and Yakovlev et al. (2001a, 2002) to interpret observations of thermal emission from isolated 
neutron stars. These authors paid special attention to the case of strong proton superfluidity 
and weak neutron superfluidity in the stellar core. Because superfluidity models have a large 
uncertainty, we consider a broader class of models without assuming from the outset that 
proton pairing is stronger than neutron pairing. In addition, attention is given to the model 
of nonstandard neutron triplet-state pairing with an anisotropic gap that vanishes along the 
quantization axis. 

2 Observational data 

The observational data on thermal emission from eleven isolated middle-aged (10 3 < t < 10 6 yr) 
neutron stars are collected in Table. In what follows, T s °° is the stellar surface temperature 
as detected by a distant observer, and t is the stellar age. The data differ from those pre- 
sented previously (see, e.g., Yakovlev et al. 2002), because they include the results of new 
observations. 

Two young objects, RX J0822-4300 and IE 1207.4-5209 (=J1210-5226), are radio-quiet 
neutron stars in supernova remnants. Two of the three oldest objects (t > 5 x 10 5 yr), 
RX J1856. 4-3754 and RX J0720. 4-3125, are also radio-quiet neutron stars. The remain- 
ing seven sources — PSR J0205+6449, the Crab pulsar (PSR B0531+21), the Vela pul- 
sar (PSR B0833-45), PSR B1706-44, PSR J0538+2817, Geminga (PSR B0633+1746), and 
PSR B1055-52 — are observed as radio pulsars. 

PSR J0205+6449 and the Crab pulsar are located in the remnants of historical super- 
novae; their ages are known exactly. 

The age of RX J0822-4300 was determined from the age of the remnant of the host 
supernova Puppis A and lies within the range t = (2 — 5) x 10 3 yr (see, e.g., Arendt et al. 
1991), with the most probable value t = 3.7 x 10 3 yr (Winkler et al. 1988). 

The age of IE 1207.4-5209 is assumed to be equal to the age of the remnant of the host 
supernova G296.5+10. According to Roger et al. (1988), this age ranges from ~ 3 x 10 3 yr 
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Table 1: Surface temperatures of isolated neutron stars 



Source 


t 


s 


Mo- 


Confi- 


References 




[10 3 yr] 


[10 6 K] 


del a > 


dence 




PSR J0205+6449 


0.82 


<1.1 


bb 


_ 


Slane et al. (2002) 


Crab 


1 


<2.0 


bb 


99.7% 


Weisskopf et al. (2004) 


RX J0822-4300 


2-5 


1.6-1.9 


H 


90% 


Zavlin et al. (1999) 


IE 1207.4-5209 


3-20 


1.4-1.9 


H 


90% 


Zavlin et al. (2003) 


Vela 


11-25 


0.65-0.71 


H 


68% 


Pavlov et al. (2001) 


PSR B 1706-44 


~ 17 


o.82±8:§i 


H 


68% 


McGowan et al. (2004) 


PSR J0538+2817 


30 ±4 


~ 0.87 


H 




Zavlin, Pavlov (2003) 


Geminga 


~ 340 


~ 0.5 


bb 


90% 


Zavlin, Pavlov (2003) 


RX J1856.4-3754 


~ 500 


<0.65 


bb 




see text 


PSR B1055-52 


~ 540 


~ 0.75 


bb 




Pavlov, Zavlin (2003) 


RX J0720.4-3125 


~ 1300 


~ 0.51 


H 




Motch et al. (2003) 



a ) Observations interpreted in terms of either a hydrogen atmosphere model (H) 
or the blackbody model (bb) 



to ~ 20 x 10 3 yr. 

The age of the Vela pulsar is assumed to lie within the range from the standard charac- 
teristic pulsar age of 1.1 x 10 4 yr to the age of 2.5 x 10 4 yr, obtained by Lyne et al. (1996) 
by analyzing the pulsar spindown with allowance for observed pulsar glitches. 

Kramer et al. (2003) estimated the age of PSR J0538+2817, t = (30 ± 4) kyr, from the 
measured proper motion of the neutron star relative to the center of the remnant of the host 
supernova S147. 

The age of RX J1856. 4-3754 was estimated by Walter (2001) from kinematic consider- 
ations and revised by Walter and Lattimer (2002). Following the latter authors, we take a 
mean value of t = 5 x 10 5 yr and choose an errorbar for t that excludes the value of t = 9 x 10 5 
yr obtained by Walter (2001). 

Zane et al. (2002) and Kaplan et al. (2002) estimated the characteristic age of RX J0720.4- 
3125 from X-ray measurements of the spindown rate of the star P. We take a mean value 
of 1.3 x 10 6 yr with an uncertainty by a factor of 2. 

The ages of the three radio pulsars — PSR B1706-44, Geminga, and PSR B1055-52 - 
are set equal to their characteristic ages with the same uncertainty factor of 2. 

For the two youngest objects (the Crab pulsar and PSR J0205+6449), only upper limits 
were placed on T s °° (Weisskopf et al. 2004; Slane et al. 2002). The surface temperatures of 
five sources — RX J0822-4300, IE 1207.4-5209, Vela, PSR B1706-44, and PSR J0538+2817 
- were determined by using neutron-star hydrogen atmosphere models (for references, see 
Table). These models give more realistic neutron star radii and hydrogen column densities 
(see, e.g., Pavlov et al. 2002) than the blackbody model. 

The pulsar PSR B0656+14 that was considered previously (see, e.g., Yakovlev et al. 2002) 
is excluded from Table. A simultaneous analysis of new X-ray and optical observations of 
the source (with the improved distance to it obtained from the parallax measurements by 
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Brisken et al. 2003) leads either to an overly small neutron-star radius (in the blackbody 
model) or to an overly small distance to the star (in the hydrogen atmosphere model); see 
Zavlin and Pavlov (2003). This makes the interpretation of thermal emission of the source 
too unreliable. 

For Geminga and PSR B1055-52, the blackbody model is more self-consistent. Therefore, 
we take the values of T s °° obtained by interpreting the observed spectra using this model. 
For PSR B 1055-52, we take T s °° from Pavlov and Zavlin (2003). 

The surface temperature of RX J1856. 4-3754 has not been determined accurately enough. 
The wide spread in T s °° obtained for different radiation models (see, e.g., Pons et al. 2002; 
Braje and Romani 2002; Burwitz et al. 2003; Pavlov and Zavlin 2003; Trumper et al. 2003) 
stems from the fact that optical and X-ray observations cannot be described by a single 
blackbody model. This may be attributed, for example, to the presence of hot spots on the 
stellar surface. Therefore, we fix only the upper limit of T s °° < 6.5 x 10 5 K that agrees with 
the value of T s °° obtained in the model of a Si-ash atmosphere (Pons et al. 2002) and in the 
model of condensed matter on the stellar surface (Burwitz et al. 2003). This limit is consistent 
also with the model of a nonuniform stellar surface temperature distribution proposed by 
Pavlov and Zavlin (2003). For the latter model, the mean stellar surface temperature is 
T s °° = 5 x 10 5 K and lies below the chosen upper limit. 

Finally, we take the surface temperature of RX J0720. 4-3125 from the paper by Motch et 
al. (2003). These authors interpreted the observed spectrum using an hydrogen atmosphere 
model of finite depth. 

For PSR J0538-4300, PSR B1055-52, and RX J0720-3125, the errors in T s °° were not 
given by the authors (see Table). In all these cases, we assume them to be equal to 20%. 

3 Models for nucleon superfluidity and neutrino emis- 
sion due to Cooper pairing of protons 

Neutron or proton superfluidity can be characterized by the critical temperature as a function 
of density, T c (p). Microscopic theories predict (see, e.g., Lombardo and Schulze 2001; for 
references, see also the review by Yakovlev et al. 1999a) the existence of singlet-state ( 1 5'o) 
neutron pairing (T cn = T cns ) in the inner crust and the outermost layers of the stellar core 
and singlet-state proton pairing (T cp ) and triplet-state ( 3 -P2) neutron pairing (T cn = T cnt ) 
in the stellar core. One should bear in mind the possibility of different components rrij of 
momentum of neutron-neutron pairs with respect to the quantization axis (\mj\=0, 1, 2) 
when triplet-state pairing is considered. A superposition of states with different mj can also 
be an energetically favored state of Cooper pairs (see, e.g., Amundsen and Ostgaard 1985; 
Baldo et al. 1992; Khodel et al. 1998, 2001). Only one type of triplet-state superfluidity 
with rrij = has commonly been assumed in calculations of neutron star cooling (except 
for the papers by Schaab et al. 1998, and Gusakov and Gnedin 2002). Below, we consider 
triplet-state pairing of neutrons with \rrij\ = and 2, because the effects of these two 
types of superfluidity on the heat capacity and the neutrino luminosity of neutron stars 
are qualitatively different. Following Yakovlev et al. (1999a), we denote the three types of 
superfluidity considered here - 1 S , 3 P 2 (^j = 0), and 3 P 2 (\mj\ =2) — by the letters A, B, 
and C, respectively. The energy gap in the neutron energy spectrum, e(p), is isotropic in 
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case A and anisotropic in cases B and C (i.e., it depends on the angle between the particle 
momentum p and the quantization axis z). In case C, the energy gap vanishes in the 
directions parallel and antiparallel to the z axis. 

Nucleon superfluidity suppresses the neutrino processes involving nucleons, changes the 
nucleon heat capacity, and triggers an additional neutrino emission mechanism related to 
Cooper pairing of nucleons (Flowers et al. 1976). The effect of neutron superfluidity C on 
the heat capacity of the matter and the neutrino reactions differs qualitatively from the 
effect of superfluidity A or B. For example, the suppression of the neutrino processes and 
the heat capacity by superfluidity C and superfluidity B or A has power-law and exponential 
character, respectively (see, e.g., Yakovlev et al. 1999a). 

Microscopic theories predict a variety of T c (p) profiles (see, e.g., Lombardo and Schulze 
2001). The peaks of T c (p) can take on values from < 10 8 K to 5 x 10 10 K. In many 
models, the peaks of T cnt (p) are lower than the peaks of T cp (p) and T cns (p), because of 
weaker nucleon-nucleon attraction in triplet-state channels. 

We use four phenomenological model profiles T c (p) of critical temperature (for both, 
neutrons and protons) in the core of a neutron star. In Fig. 1, these models are denoted by a, 
b, c, and d. The chosen T c (p) profiles are similar and differ only in height (maximum value): 
T max = 10 io } 4.0 x 10 9 , 8.0 x 10 8 , and 8.0 x 10 7 K (models a, b, c, and d). Superfluidities a, 
b, c, and d will be called strong, moderately strong, moderate, and weak, respectively. The 
chosen models are consistent with theoretical calculations of T c (p). The T c (p) curves have 
steep slopes at p > p^, where p^ is the threshold density at which the direct Urea process 
is open (see below). 

Below, we will denote combined nucleon superfluidity by a/5, where a is one of neutron 
triplet-state (type B or C) superfluidity models (a, b, c, or d), and (3 is one of proton 
singlet-state superfluidity models (a, b, c, or d). 

Note a large uncertainty in the neutrino emissivity Q p due to Cooper pairing of protons. 
In the non-relativistic approximation (Yakovlev et al. 1999b), Q p oc ( p , where £ p = Cy p , 
and c Vp ~ 0.08 is the vector constant of neutral current of protons; the latter constant is 
numerically small and leads to unreasonably low values of Q p . For comparison, the neutrino 
emissivity Q Q due to triplet-state pairing of neutrons is proportional to ( n = c\ n + 2c An = 
4.17, where cy n = 1 and CA n = —1.26 are the vector and axial vector constants of neutral 
neutron current, respectively. According to Kaminker et al. (1999), the inclusion of the 
relativistic correction that contains the axial vector constant ca p = 1.26 of neutron current 
of protons can greatly (by a factor of 10 to 50) increase the constant ( p (and the emissivity 
Q p ) compared to the non-relativistic value of ( p = Cy p = 0.0064. 

We used this value of £ p , enhanced by relativistic effects, in our previous simulations of 
neutron star cooling. On the other hand, while studying the cooling of stars with density 
dependent critical temperatures of protons T cp (p) (see, e.g., Kaminker et al. 2002), we re- 
stricted ourselves to the models for strong proton superfluidity (similar to model a). Such 
a superfluidity arises at early cooling stages. At these stages, the neutrino emission due to 
proton pairing cannot compete with other neutrino processes and plays no special role. In 
a cooler star, this neutrino emission is generated only in a small volume and weakly affects 
the cooling. Thus, in the cooling scenarios with strong proton superfluidity, considered pre- 
viously, the emission due to Cooper pairing of protons was unimportant (as well as the exact 
value of ( p ). 
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Figure 1: Model density profiles of the critical neutron and/or proton temperatures in the core of a 
neutron star. The vertical dot-dashed line indicates the threshold density at which the direct Urea process 
becomes operative. 



In this paper, we consider (among others) models of moderate proton superfluidity in 
which the emission due to proton pairing can appreciably affect the cooling, so that the 
value of ( p is important. As noted, for instance, by Yakovlev et al. (1999b) and Kaminker 
et al. (1999), the constant ( p can De affected not only by relativistic effects, but also by 
renormalization due to many-body effects in nucleon matter. This renormalization for the 
process in question has not yet been performed. Carter and Prakash (2002) gave an example 
of a similar renormalization of the constant of the axial vector current. For the sake of 
definiteness, we perform calculations by choosing the renormalized value of £ p — 1- The 
sensitivity of our calculations to the value of ( p is described in Section 5.3 (also see Fig. 7). 



4 Cooling of stars with strong proton superfluidity 

Let us compare the observational data with our calculated cooling curves (T s °°(t) profiles). 
The calculations were performed using a code described by Gnedin et al. (2001). As in 
previous papers (mentioned in the introduction), we consider the models of neutron stars 
whose cores are composed of neutrons n, protons p, and electrons e. We use a moderately 
stiff equation of state in the stellar core proposed by Prakash et al. (1988) (model I with 
the compression modulus of symmetric nucleon matter at saturation K = 240 MeV). The 
maximum mass of a stable neutron star, for the chosen equation of state, is M = 1.977M Q 
(at a radius of R = 10.754 km and a central density of p c = 2.575 x 10 15 g cm -3 ). This 
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equation of state opens the intense direct Urea process of neutrino emission (Lattimer et al. 
1991) at densities p above the threshold density pn = 7.851 x 10 14 g cm -3 , i.e., in stars with 
masses M > M D = 1.358M Q . The radius of a star with the threshold mass M D is R = 12.98 
km. 

The thermal evolution of a neutron star consists of three stages: 

(1) the stage of thermal relaxation of the inner stellar layers (t < 100 yr); 

(2) the subsequent stage of neutrino cooling (10 2 < t < 10 5 yr) of a star with an 
isothermal core via neutrino emission from the inner layers (mainly from the core); 

(3) the final stage of photon cooling (t > 10 5 yr) via photon emission from the stellar 
surface. 

The cooling theory for non-superfluid stars cannot explain the entire set of observational 
data (see, e.g., Kaminker et al. 2002). However, the theory can be reconciled with the 
observations by taking into account nucleon superfluidity. According to Kaminker et al. 
(2001), it is sufficient to assume the presence of strong proton superfluidity and weak neutron 
superfluidity in stellar cores. 

Figure 2 shows cooling curves for neutron stars of different masses with weak neutron 
superfluidity d and strong proton superfluidity a. Such weak neutron superfluidity switches 
on only at the photon cooling stage. Therefore, the type of weak neutron superfluidity (B 
or C) does not affect the cooling of middle-aged stars. The family of cooling curves for stars 
with masses M > M fill in the hatched region. All of the observed sources fall within this 
region; i.e., they can be interpreted in terms of the assumed superfluidity model. 

As shown Kaminker et al. (2002), strong proton superfluidity (with weak neutron super- 
fluidity or with normal neutrons) gives rise to three types of cooling neutron stars. 

Low- mass stars cool down very slowly (more slowly than low-mass non-superfluid stars). 
Cooling curves for such stars weakly depend on their mass, the equation of state in their 
cores, and proton superfluidity model (on the specific form of the T cp (p) profile as long as 
superfluidity in the stellar core is strong enough, T cp (p) > 4 x 10 9 K). The upper boundary 
of the hatched region in Fig. 2 is the cooling curve for a star with a mass of M — 1.35M ; 
it is almost indistinguishable from the cooling curve for a star with M = 1.1 M Q and agrees 
with the observations of four sources, RX J0822-4300, IE 1207.4-5209, PSR B1055-52, 
and RX J0720. 4-3125, the hottest ones for their ages. These sources will be considered as 
low-mass neutron stars. 

High-mass neutron stars cool down very rapidly via intense neutrino emission generated 
by the direct Urea process in the inner stellar core. At high densities (p > 10 15 g cm -3 ), 
proton superfluidity weakens (Fig. 1) and ceases to suppress the neutrino emission. The 
cooling curves for such stars weakly depend on their mass, the equation of state, and the 
proton superfluidity model. They almost coincide with the cooling curves for high-mass 
non-superfluid stars. All of the observed isolated neutron stars are much hotter than the 
stars of this type. 

Finally, medium-mass stars cool down moderately rapidly. Their cooling strongly de- 
pends on the mass, the equation of state, and proton superfluidity model. By varying the 
stellar mass, we can obtain a family of cooling curves that fill the space between the cooling 
curves for low-mass and high-mass stars. We consider the sources PSR J0205+64, Vela, 
PSR B1706-44, PSR J05438+2817, Geminga, and RX J1856.4-3754 as medium-mass neu- 
tron stars. 
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2 4 6 



lg t (yrs) 

Figure 2: Comparison of the observations (see Table) with the cooling curves for neutron stars with 
masses from 1.1 to 1.977-M© (indicated near the curves) for weak neutron superfluidity d and strong proton 
superfluidity a. The region filled by cooling curves for stars of different masses is hatched. 



5 Cooling of neutron stars with combined nucleon su- 
perfluidity 

Figures 3-6 show cooling curves for neutron stars with different superfluidities of neutrons 
a and protons /3 (a, (3 = a, b, c, or d). Neutron superfluidity is of type B. We consider all 
the possible combinations of neutron and proton superfluidities. For each combination a/3, 
we show the upper cooling curve of a low- mass star (M = 1.1M , with the central density 
p c = 6.23 x 10 14 g cm~3) and the lower cooling curve of a high-mass star (M = M max ). 
The lower curve is virtually independent of the models of superfluidity a(3 (see the previous 
section). The region between the upper and lower curves (similar to the hatched region 
in Fig. 2) can be filled by cooling curves of medium-mass stars and is accessible by an 
assumed model of superfluidity a(3. As in Fig. 2, we show also the observational data. 
Nucleon superfluidity models can be constrained by comparing accessible T s °° regions with 
the observational data. 

Each of Figs. 3-6 consists of two panels: In panel (a) neutron superfluidity a is fixed, and 
cooling curves are given for all four proton superfluidity models; in panel (b) proton super- 
fluidity (3 is fixed, and cooling curves are given for all four neutron superfluidity models. By 
comparing panels (a) and (b), we can trace the change in cooling when proton superfluidity 
is replaced by neutron superfluidity (and vice versa). 
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Figure 3: (a) Cooling of stars of two masses, M = 1.977M and M = 1.1M Q , for neutron 
superfluidity a and different proton superfluidity models (a, b, c, or d); (b) the same for 
proton superfluidity model a and different neutron superfluidity models [a, b, c, or d). 
Neutron superfluidity B is employed in all the cases. Theoretical curves are compared with 
the observations. The cooling of a star with M = 1.977-M© does not depend on superfluidity 
model. 

5.1 Fixed proton superfluidity 

Let us choose a proton superfluidity model (f3 = a, b, c, or d) and consider the dependence 
of the upper cooling curves on the models of neutron superfluidity a in Figs. 3-6. 

The cooling curves b(3 run below the curves a(3 because of the neutrino emission due to 
Cooper pairing of neutrons. This emission is significantly suppressed in the models with 
strong neutron superfluidity a (see, e.g., Yakovlev et al. 1999a, 1999b, 2001b). All other 
neutrino reactions involving neutrons and the neutron heat capacity are fully suppressed by 
neutron superfluidity a or b. The difference between the cooling curves a(3 and b(3 depends 
on the model for proton superfluidity (3. Thus, for example, as we go from the model f3 = a 
to (3 = b and then to the model of moderate proton superfluidity (3 = c, the contribution 
of the neutrino emission due to Cooper pairing of protons to the neutrino luminosity of 
the star increases (and becomes dominant for (3 = c). Indeed, the neutrino emission due 
to Cooper pairing affects most strongly the cooling at moderate critical temperatures of 
nucleons, T c ~ 2 x (10 8 - 10 9 ) K (see, e.g., Yakovlev et al. 1999b, 2001b). As a result, the 
difference between the cooling curves a(3 and b(3 in Figs. 3-5 steadily decreases as we go from 
the model (3 = a to (3 = b and (3 = c. At the same time, the accessible theoretical cooling 
regions disagree with the observations more and more. 

For weak superfluidity (3 = d (Fig. 6), protons remain normal as long as t < 10 5 yr, 
until the onset of the photon cooling stage. In this case, the main neutrino process involving 
protons is the bremsstrahlung in proton-proton collisions. Since the contribution of this 
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Figure 4: Same as in Fig. 3 but for fixed superfluidity model b of neutrons (a) or protons 
(b). 

process to the neutrino emission is much smaller than the contribution of Cooper pairing of 
protons in model c, the difference between the curves ad and bd again increases. As in Fig. 
3 (curves aa and ba), it is mainly determined by the more intense neutrino generation due 
to Cooper pairing of neutrons in model b than in model a. As a result, the accessible region 
of stellar surface temperatures for combined superfluidity ad (as for superfluidity da, cf. Fig. 
6 with Figs. 2 and 3) agrees with the observations. 

The curves af3 and bf3 approach one another at the photon cooling stage (t > 10 5 yr). In 
this case, the influence of neutron superfluidity a or b on the cooling of the star manifests 
itself mainly in strong suppression of the neutron heat capacity. As a result, the heat capacity 
of the star is determined by the heat capacity of protons (also suppressed by superfluidity 
(3) and electrons. 

For superfluidity models c/3, the neutrino emission due to Cooper pairing of neutrons 
is especially efficient. Accordingly, the cooling curves c(3 in Figs. 3-6 run well below the 
curves a(3 and bf3 and do not differ too much from one another. In particular, all of the 
curves c(3 describe the sharp speedup in cooling at t ~ 300 yr associated with the switch- 
on of the neutrino emission due to neutron pairing. We can see that the accessible surface 
temperatures obtained for the superfluidity models c(3 lie well below most of the observational 
data points. 

For the model of weak neutron superfluidity d, the d(3 curves almost coincide with the 
cooling curves for normal neutrons. Differences arise only at the photon cooling stage 
(t > 10 5 yr) from the partial suppression of the neutron heat capacity. However, at the 
neutrino cooling stage, neutron superfluidity d has not yet set in. Therefore, all of the d(3 
cooling curves lie above the c(3 curves. At t > 10 5 — 10 6 yr, the dj3 cooling curves for any (3 
run above the curves a(3 and b(3 due to the strong suppression of the neutron heat capacity 
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by superfluidities a and b. Finally, the cooling curve for the model of combined superfluidity 
dd is close to the standard cooling curve for non-superfluid low-mass (M < Md) neutron 
stars. This cooling curve disagrees with the observations of many neutron stars (both the 
hottest and coolest ones for their ages). 




2 4 6 2 4 
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Figure 5: Same as in Fig. 3 for superfluidity model c of neutrons (a) or protons (b). 



5.2 Fixed Neutron Superfluidity 

Let us choose a neutron superfluidity model (a = a, b, c, or d) and consider the dependence 
of the upper cooling curves on proton superfluidity models (3. A comparison of panels (a) 
and (b) in Figs. 3-6 reveals a qualitative similarity between cooling curves of low-mass 
(M = 1.1M ) stars where neutron and proton superfluidities are inverted (i.e., for models 

a/3 and (3a). 

The quantitative differences between cooling curves a(3 and [3a are attributable to dif- 
ferent neutron and proton number densities in the cores of neutron stars and to different 
types of neutron (triplet-state) and proton (singlet-state) pairing. This results in a slightly 
asymmetric effect of neutrons and protons on the neutrino luminosity and the heat capac- 
ity (see, e.g., Yakovlev et al. 1999a). Thus, for example, at temperatures T slightly below 
T c , the neutrino emissivity due to neutron pairing is approximately an order of magnitude 
higher than that due to proton pairing (even for the chosen constant ( p = 1). Therefore, 
the cooling curves b((3 = a, d) and c(f3 = a, b, d) (panel (a) in Figs. 4 and 5) lie below the 
"inverted" cooling curves (a = a,d)b and (a = a,b,d)c (panel (b) in the same figures). 
On the other hand, the curve ad (Fig. 3a) at the photon cooling stage (t > 10 5 yr) runs 
below the curve da (Fig. 3b). This is because the heat capacity of the neutron star core is 
stronger suppressed by neutron superfluidity a = a than by proton superfluidity (3 = a. In 
other cases, the inversion of neutron and proton superfluidities leads to qualitatively similar 
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Figure 6: Same as in Fig. 3 for superfluidity model d of neutrons (a) or protons (b). 

(roughly symmetric) cooling curves in Figs. 3-6. For high-mass neutron stars (M > Md), 
this symmetry was found by Levenfish et al. (1999) in their simplified cooling calculations 
for stars with constant critical temperatures of neutrons and protons over the stellar core. 

A comparison of the upper cooling curves with the observations in Figs. 3-6 shows that 
there are only two models of combined nucleon superfluidity that are consistent with the set 
of observational data. They include the model da discussed in Section 4 and the "inverted" 
model ad (Figs. 3 and 6). In other words, one (neutron or proton) superfluidity must be 
weak, while the other must be strong. Other models are unable to simultaneously explain 
the observational data, primarily for four neutron stars (RX J0822-4300, IE 1207.4-5209, 
PSR B1055-52, and RX J0720.4-3125), the hottest ones for their ages. 

Varying nucleon superfluidity models, we can constrain critical temperatures of nucleons 
at which the theory agrees with the observations. In general, the following conditions must be 
satisfied simultaneously: either T™ ax < 2 x 10 8 K and T™ ax > 4 x 10 9 K, or T c ™ ax > 5 x 10 9 K 
and T™ ax < 2 x 10 s K. 

The models of moderate neutron and/or proton superfluidity in the cores of neutron 
stars with peak temperatures T™ ax and/or T™ ax in the range ~ (2 x 10 8 — 4 x 10 9 ) K are 
inconsistent with the observations of neutron stars hottest for their ages. We can show 
that this conclusion is valid for a much broader class of nucleon superfluidity models than 
those used here. Nevertheless, there is a narrow region of nucleon superfluidity parameters 
at which the combination of strong nucleon superfluidity of one type and moderate nucleon 
superfluidity of another type can be reconciled with the observations (for details, see Gusakov 
et al. 2004). 
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Figure 7: Cooling curves of a low-mass star (M = 1.1M ) for three models of nucleon 
superfluidity (db - long dashes, dc - dash-and-dot curves, be - short dashes) in comparison 
with the observations. Lower and upper curves for each superfluidity model are computed, 
respectively, with the renormalized (C p = 1) and non-renormalized constant in the expression 
for neutrino emissivity due to Cooper pairing of protons. 

5.3 On the constant in the expression for the neutrino emissivity 
due to Cooper pairing of protons 

Let us briefly discuss the sensitivity of cooling curves to the constant ( p in the expression 
for the neutrino emissivity due to proton pairing. Recall that the value of ( p that includes 
many-body effects is known poorly. In our calculations, we used the (renormalized) value 
C„ I- 

For example, Fig. 7 shows cooling curves of a low-mass neutron star for three models of 
neutron and proton superfluidity (db, dc, and be). As shown above, the neutrino emission due 
to Cooper pairing is especially important in low-mass stars. As everywhere in this section, we 
consider neutron superfluidity of type B. The lower of two curves for each superfluidity model 
is calculated with the renormalized constant ( p = 1, while the upper curve is calculated with 
the non-renormalized constant (but taking into account relativistic effects; see Kaminker et 
al. 1999). 

In model db, proton superfluidity b is moderately strong and appears at an early cooling 
stage. The neutrino emission due to proton pairing plays a relatively minor role, and the 
exact value of Q p weakly affects the cooling. 

In model dc and, especially, in model be, moderate proton pairing f3 = c results in intense 
neutrino emission and appreciably speeds up the cooling. In these cases, the cooling curves 
are most sensitive to £ p . However, as seen from Fig. 7, the employed variations of £ p cannot 
lead to agreement of cooling curves dc and be with the observations and, hence, do not affect 
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our conclusions. We believe the renormalized value ( p = 1 to be more realistic than the non- 
renormalized value. The existing uncertainty in ( p introduces an uncertainty in the cooling 
theory. In particular, for the non-renormalized value of ( p , the approximate symmetry of 
cooling curves relative to the inversion of nucleon superfluidity models {af3 # (3a) noted 
above is much less pronounced than for the renormalized value (see also Yakovlev et al. 
1999a). It is possible that the choice of ( p may be important in the future for reconciling 
theory with observations. 

6 Two types of triplet-state pairing of neutrons 

Let us compare the effect of the two types of neutron triplet-state superfluidity (B and C) on 
cooling of neutron stars. Clearly, significant differences might be expected for strong neutron 
superfluidity. As follows from the results of the previous section, strong neutron superfluidity 
of type B (model a) and weak proton superfluidity (model d) can ensure agreement between 
the theory and the observations. Let us consider this case in more detail. Figure 8 shows 
the cooling curves for neutron stars of different masses with superfluidity ad: In panels (a) 
and (b), we take neutron pairing of type B and C, respectively. 

According to Fig. 8a (as for model da in Fig. 2), we can identify the same three types of 
cooling neutron stars: low-mass, slowly cooling stars; high-mass, rapidly cooling stars; and 
medium-mass stars with a moderate cooling rate. 

The cooling of neutron stars with neutron superfluidity of type C was first calculated 
by Schaab et al. (1998). However, these authors used an oversimplified description of the 
effects of superfluidity on neutrino reactions. More accurate calculations were performed 
by Gusakov and Gnedin (2002), who compared the results obtained for superfluidities of 
types B and C. The authors used the approximation of critical temperatures T cp and T cnt 
constant over the stellar core. Calculations indicate that, in many cases, cooling curves do 
not change if the actual T cp (p) and T cnt (p) profiles are replaced by effective constant critical 
temperatures close to T cp (p c ) and T cnt (p c ) at the stellar center (p = p c ). This approximation 
is valid if T c (p) is a smooth function of p near the stellar center (e.g., in low-mass stars). 

Gusakov and Gnedin (2002) showed that neutron superfluidity of type C speeds up the 
neutron star cooling (compared to superfluidity of type B). This is caused by the power-law 
suppression of the neutrino emissivity by superfluidity C (in contrast to the exponential 
suppression in case B; see, e.g., Yakovlev et al. 1999a; Gusakov 2002). Our calculations (Fig. 
8) indicate that the above conclusion remains valid in a more realistic approach, taking into 
account variations of the critical temperatures T cnt (p) and T cp (p) over the stellar core. 

Cooling curves for low- and medium-mass stars in Fig. 8b lie well below cooling curves 
for stars of the same masses in Fig. 8a. On the other hand, cooling curves for high-mass 
stars (M > 1.55M Q ) in both panels of the figure almost coincide for the obvious reason: the 
critical temperatures are low in the central regions of these stars (T cnt < 10 8 K, see Fig. 1), 
so that superfluidity ceases to affect the cooling. 

Thus, according to Fig. 8, strong neutron superfluidity of type C disagrees with the 
observations of the neutron stars hottest for their ages. This superfluidity is still too weak to 
completely suppress the modified Urea process in a low-mass star, thereby making the star 
hotter. Of course, the theory can be reconciled with the observations by choosing the model 
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of stronger superfluidity C. Our calculations indicate that this requires a T cnt (p) profile with 
the maximum T™f x ~ 10 11 K. However, such a strong triplet-state pairing seems unrealistic. 
Khodel et al. (1998, 2001) gave theoretical arguments against the appearance of superfluidity 
of type C in neutron stars. 
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Figure 8: Same as in Fig. 2 for model a of strong neutron superfluidity of type B (a) or C 
(b) and model d of weak proton superfluidity (of type A). 



7 Conclusions 

We have analyzed the cooling of neutron stars with triplet-state pairing of neutrons and 
singlet-state pairing of protons in stellar cores. Our results are summarized below. 

(1) Cooling curves are qualitatively symmetric relative to the inversion of neutron and 
proton superfluidity models. For low-mass stars (M < Md), this symmetry is largely de- 
termined by the dimensionless constant ( p ~ 1 in the expression for the neutrino emissivity 
due to Cooper pairing of protons. At ( p -C 1, obtained without the renormalization of ( p 
with account for many-body effects, the symmetry is much less pronounced than at ( p ~ 1 
expected after the renormalization. 

(2) Two types of neutron and proton superfluidity models are consistent with the ob- 
servations of thermal emission from isolated neutron stars. First, they are the models of 
strong (type A) proton superfluidity and weak (or absent) (type B or C) neutron super- 
fluidity with maximum (over the stellar core) critical temperatures T™ ax > 4 x 10 9 K and 
T™f x < 2 x 10 s K. Second, they are the models of strong (type B) neutron superfluidity 
and weak (or absent) proton superfluidity with T™ ax > 5 x 10 9 K and T™ ax < 2 x 10 8 K. 
The models of the first type seem more realistic. Note, in particular, a recent paper by 
Schwenk and Friman (2004) who predicted a weakening of triplet-state pairing of neutrons 
by many-body effects. 
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(3) Models of moderate (type B or C) neutron superfluidity and/or moderate (type A) 
proton superfluidity with maximum critical temperatures T™f x and T™ ax in the range from 
~ 2 x 10 8 to ~ 4 x 10 9 K are inconsistent with the observations of primarily the hottest 
young neutron stars RX J0822-4300 and IE 1207.4-5209. However, the agreement between 
the cooling theory and the observations is also possible for a narrow range of parameters 
of moderate neutron superfluidity (T™ x ~ 6 x 10 8 K). This possibility was considered by 
Gusakov et al. (2004). 

(4) Strong triplet-state neutron pairing of type C can appreciably speed up the cooling 
of middle-aged (10 3 < t < 10 5 yr) neutron stars compared to pairing of type B for the 
same T cn t(p) profiles. For strong neutron pairing of type C {T™ x > 5 x 10 9 K) and weak 
proton pairing (T™ ax < 2 x 10 8 K), the theory cannot be reconciled with the observations of 
RX J0822-4300, IE 1207.4-5209, PSR B1055-52, and RX J0720.4-3125, the hottest sources 
for their ages. For any models of strong or moderate neutron pairing, the transition from 
pairing of type B to pairing of type C just enhances the difference between the theory and 
the observations. 

Our analysis is simplified, because we have considered only the cores of neutron stars com- 
posed of neutrons, protons, and electrons (disregarding the possible existence of hyperons, 
pion and kaon condensates, or quark matter). Moreover, we have chosen only one equation 
of state in the cores of neutron stars and similar profiles of the critical temperatures T c (p) 
of neutrons and protons in the stellar cores. Varying the equation of state (for the matter 
composed of neutrons, protons and electrons) leads to shifting the threshold of the direct 
Urea process (to changing p D and M D ). Varying T c (p) profiles (but retaining their general 
shape) at high T c max > 2 x 10 9 K leads to shifting the characteristic values of p at which 
superfluidity weakens and ceases to suppress the intense neutrino emission. Both effects 
shift boundary masses that separate the three types of cooling neutron stars (Kaminker et 
al. 2002), but do not alter our main conclusions. It is significant that the simplest model 
of neutron stars with strong proton superfluidity (even without neutron superfluidity) is 
capable of explaining the available observations. 

It should be noted that the cooling of neutron stars also depends on (singlet-state) neutron 
pairing in the inner stellar crust, on the magnetic field in the outermost stellar layers, and 
on the presence or absence of a surface layer of light elements (see, e.g., Potekhin et al. 
2003). In general, however, these factors have a weaker effect on the cooling than nucleon 
superfluidity in the cores of neutron stars considered here. We have disregarded them by 
restricting our analysis to superfluidity effects in the stellar cores. Our cooling code allows us 
to include these factors, and it can be necessary for interpreting individual sources, primarily 
the objects hottest for their ages (see, e.g., Potekhin et al. 2003). The cooling of neutron 
stars can also depend on internal reheating mechanisms related, for example, to the viscous 
dissipation of differential rotation (see, e.g., Page 1998a, 1998b). We emphasize that these 
reheating mechanisms are model dependent. On the other hand, the available observations 
can be interpreted without invoking them. 

Note that the surface temperatures of neutron stars T s °° are difficult to determine from 
observational data (see, e.g., Pavlov et al. 2002). For this purpose, one needs reliable obser- 
vational data and theoretical models of neutron star atmospheres. The existing values of T s °° 
can change appreciably, which can strongly affect the interpretation of the observations, es- 
pecially, of RX J0822-4300, IE 1207.4-5209, PSR B1055-52, and RX J0720.4-3125. Future 
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observations of thermal emission from isolated neutron stars will be crucial for understanding 
superfluid properties of dense matter in stellar cores. 
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